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ABSTRACT 


The  purpose  of  this  note  is  to  examine  the  relationships 
of  the  Kalman  filters  which  can  be  used  to  estimate  the 
linearized  state  of  two  different  state -variable  repre¬ 
sentations  of  the  same  nonlinear  system.  It  is  shown 
that,  to  within  first  order  effects,  the  choice  of  the  co¬ 
ordinate  system  for  the  construction  of  the  Kalman  filter 
is  immaterial. 
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1.  INTRODUCTION 


It  is  often  desirable  to  estimate  the  state  of  a  re-entry  vehicle,  satellite, 
etc.  by  observing  certain  output  signals,  using  a  radar,  in  the  presence  of 
noise.  In  general,  the  equations  of  motion  of  the  system  are  nonlinear;  the 
nonlinear  terms  arise  due  to  the  inverse  -  square  law  and  the  nonlinear  de¬ 
pendence  of  the  drag  force  on  the  vehicle  velocity. 

In  such  problems,  the  fact  that  one  is  using  a  radar  to  make  observations 
fixes  the  coordinate  system  in  which  the  output  signals  are  measured.  Thus, 
one  usually  observes 

a)  range, 

b)  azimuth,  and 

c)  elevation 

in  the  presence  of  measurement  noise,  which  is  in  general  assumed  to  be 
white. 

The  prediction  of  the  future  system  trajectory  necessitates  the  estima¬ 
tion  of  the  remaining  state  variables  of  the  system.  The  estimation  of  these 
variables  necessitates  the  definition  of  a  coordinate  system.  For  example, 
for  a  radar-based  coordinate  system  we  must  estimate 


a) 

range 

d) 

range  -  rate 

b) 

azimuth 

e) 

azimuth -rate 

c) 

elevation 

f) 

elevation  -  rate 

On  the  other  hand,  it  may  be  desirable  to  estimate  the  state  variables 
using  a  different  coordinate  system.  For  example,  one  may  wish  to  use  an 
inertial  coordinate  system  or  a  vehicle -based  coordinate  system.  The 
choice  of  the  coordinate  system  fixes  the  (nonlinear)  differential  equations 
of  the  chosen  state  variables  (equations  of  motion)  and  the  transformations 
linking  the  state  variables  to  the  output  variables  (range,  azimuth,  and 
elevation). 

Since  we  are  dealing  with  one  and  the  same  physical  system,  the  various 
equations  are  related.  In  point  of  fact,  given  any  one  state  variable  repre¬ 
sentation  one  can 

a)  obtain,  via  dynamic  linearization  techniques,  a  set  of 
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linear  perturbation  differential  equations,  and 
b)  construct  the  corresponding  Kalman  filter  for  the 
estimation  of  the  perturbation  variables. 

It  follows  that  each  coordinate  system  yields  a  different  linearized  Kalman 
filter.  The  natural  question  that  arises  is:  is  there  any  advantage  in  choos¬ 
ing  any  one  particular  state  variable  representation? 

It  is  shown  in  this  report  that  if  the  linearized  equations  adequately  rep¬ 
resent  the  system  behavior,  then  any  choice  of  a  coordinate  system  yields 
the  same  answer.  To  accomplish  this  we  proceed  as  follows.  In  section  2 
we  define  two  different  state  variable  representations  of  the  same  system 
and  with  the  same  output  variables.  In  section  3  we  state  the  mathematical 
relations  which  link  the  two  state  variable  representations.  In  section  4 
we  obtain  the  linearized  perturbation  equations  by  linearizing  about  a 
nominal  trajectory.  In  section  5  we  derive  the  equations  of  the  Kalman 
filters  which  can  be  used  to  estimate  the  perturbations  of  the  state  from  the 
nominal  trajectory;  we  then  show  that  it  is  immaterial,  to  within  first-order 
terms,  which  state  variable  representation  is  used  to  construct  the  Kalman 
filter. 
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2.  STATEMENT  OF  THE  PROBLEM 


Let  us  suppose  that  we  are  given  a  nonlinear,  time -invariant,  dynamical 
system  Sj  with  state  vector  x  and  output  vector  related  by  the  equations 


S 


1 


k  = 

1  =  hj(x) 


(2.  1) 
(2.2) 


We  suppose  that  this  system  is  completely  observable. 

Now  let  us  suppose  that  we  represent  the  same  system  by  using  a  differ¬ 
ent  coordinate  system  in  which  the  state  vector  is  _z  and  the  output  vector 
is  still  We  suppose  that  the  vectors  _z  and  x  are  related  by 

z  =  £(x)  (2.3) 

such  that  g(.  )  is  differentiable,  one-to-one,  and  onto  for  all  x.  Thus,  we 
suppose  that  for  every  x  there  is  a  unique  z  and  vice-versa.  In  other 
words,  the  inverse  transformation  (.  )  exists  and 


-1,  V 

x  =  £  (_z) 


(2.4) 


is  well  defined,  for  all  x  and  _z. 

In  the  new  coordinate  system  the  dynamical  system  is  described  by  the 
relations 


f  i  =  12(Z) 

{  I  =  hZ(iO 


(2.  5) 
(2.  6) 


We  next  suppose  that  we  cannot  observe  the  output  vector  ^y(t)  in  the 
absence  of  noise.  So  we  let  v(t)  denote  a  vector -valued  white  noise  process, 
with  zero  mean  and  covariance  matrix  R(t),  i.  e. 

E  |  v(t)j  =  0  (2.7) 

COV  £v(t);  v(t)]  =  E  £  v(t) v* (t )j-  =  R(t)  6  (t  -  t)  (2.8) 


The  basic  problem  is  to  estimate  the  state  x(t)  (or  _z(t)  )  by  observing 
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the  signal 


over 


w(t)  =  x(T)  +  Z<T) 
period  of  time 

t  ^  T  <  t 
O 


(2.  9) 


(2.  10) 
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3.  SOME  MATHEMATICAL  RELATIONS 


In  this  section  we  shall  state  some  self-evident  mathematical  relations 
that  relate  the  two  system  representations  S^  and  S^. 

First  we  start  with  the  algebraic  relation  (2.  3),  i.  e. 


Z  =  £.(x) 


(3.  1) 


We  differentiate  with  respect  to  time  and  use  the  chain  rule  to  obtain 

.  I  I  .  '  dM  l 

£  =j-5x- 1?  =  —  l (— > 

where  (  d j*/  5x)  is  the  Jacobian  matrix  of  g(.  ).  Now  since 

A  =  lz(i) 

then,  in  view  of  Eq.  (3.  1),  we  have 

i  =  i2(g(*) ) 

From  Eqs.  (3.2)  and  (3.4)  we  deduce  the  relation 


(3.2) 


(3.3) 


(3.4) 


I  \ 

=i2 (£(*))  all  x  (3.5) 

Similarly,  from  Eqs  (2.  2)  and  (2.  6)  we  obtain,  in  view  of  Eq.  (3.  1), 
the  relation 


hj  (x)  =  h2(^(x))  for  all  x 
Let  us  now  define  the  matrix  G  (x)  by 


(3.  6) 


G(x)=  |-a- 


and  the  matrix  £(z)  (the  Jacobian  matrix  of 

r  I  h&~1\ 

!<*>  =*  TS 


(3.  7) 


(z))by 


(3.8) 
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Since 


z  =  G(x)x 


(3.9) 


and 


x  =  r(z)z  (3.  10) 

it  follows  that 

G  (x)  r  (z)  =  1^  for  all  x  and  _z  .  (3.  1  1) 

where  I  is  the  identity  matrix. 
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4.  DYNAMIC  LINEARIZATION 


Let  t  be  some  initial  time.  Let 
o 

x*(tQ)  (4.  1) 

denote  some  nominal  initial  state  for  the  system  Sj.  Let 

z*(t0)  =  £(x*(to)  )  (4.2) 

denote  the  nominal  initial  state  for  the  system  S^.  Let 

X*  (t)  ;  t  >  to  (4.3) 

denote  the  nominal  trajectory  (solution)  of  the  differential  equation  (2.  1) 
and  let  z  (t)  ,  given  by, 

z'  (t)  =  _g(  x'  (t)  )  (4.  4) 

denote  the  nominal  trajectory  (solution)  of  the  system  S^. 

We  shall  denote  by  _?  (t)  small  perturbations  about  the  trajectory 
x  (t)  so  that  the  actual  state  of  x(t)jis  given  by 

x(t)  =  x'  (t)  +  J  (t)  (4.  5) 

We  shall  denote  by  J(t)  small  perturbations  about  the  trajectory  z  (t) 
so  that  the  actual  state  of  S2^£(t)^is  given  by 

z(t)  =  z  (t)  +  J  (t)  (4.  6) 


Then  for  the  system  Sj  we  have 


x*(t)  +  |(t)  =I1(x*(t))  + 


-1  I 

^x~|x*(t) 


5(t)  +0(5) 


(4.7) 


where  o(§)  are  the  (vector-valued)  higher  order  terms, 
matrix  Fj(x) 


We  define  the 


(4.  8) 
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and  the  matrix 


Fj  (t)  =  Fj(x  (t))  = 


*1 


dx  x  (t) 


(4.  9) 


By  neglecting  the  higher  order  terms  in  Eq.  (4.  7)  we  have 


1  =  F*(t)| 


(4.  10) 


In  a  similar  manner  we  have  for  the  system 


i"'(t)  +  j(t)  =  f2(z"(t))  + 


ai2 


Ti  /- 


s*<t>  + 


(4.  11) 


Define  the  matrices 


F2(z)  2 


siz 


and 


i  —2  (— *  <t)  ) 


(4.  12) 


(4.  13) 


By  neglecting  the  higher  order  terms  o(J)  in  Eq.  (4.  11)  we  have 


j  ~  f2'  (t)  J 


(4.  14) 


We  let  y  (t)  denote  the  nominal  output  of  the  system.  Let  J]  (t)  denote 
the  perturbations  about  the  nominal  output  so  that  the  actual  output  y(t), 
for  both  systems,  is 


x(t)  =  X  (t)  +  2(t) 


(4.  15) 


From  Eqs.  (2.  1)  and  (4.  15)  we  have 


X'  (t)  +  £(t)  =  2i !  (iS  (t)  )  + 


S  hj 


x  / 


x*(t)  JW  + °i(l> 


(4.  16) 
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while  from  Eqs.  (2.  6)  and  (4.  15)  we  have 


+  T](t)  =  h2(z»)  + 


sVi 

7T-]z*(t)I(t)  +  £lU) 


Define  the  matrices: 


Hj  ( x  )  § 


9hj 


1  w 


t  Hj ( x  (t)  )  = 


Shi' 


3  x  j  x'(t) 


h2(£)  i 


Sh2 


)  z  / 


H2  (t)  §  H2  (  z  (t)  )  2 


z"(t) 


Clearly 

J2  (t)  =  H*  (t)  |(t)  +  oj  (£)  =  H*  (t)  J(t)  +  o^J) 

If  we  neglect  the  higher  order  terms  ^  (§)  and  ^(J)  we  have 
n  (t)  =  H;(t)i(t) 
n(t)  =  H2(t)  j(t) 

Recall  that  (see  Eq.  (2.  29)  )  the  observed  signal  is 
w(t)  =  y(T)  +  v  (t) 

so  that 

w(t)  =  %  (t)  +  ti(t)  +  v(t) 

j/ 

Since  y  (t)  is  a  deterministic  and  computable  signal,  we  can  think 
the  observed  signal  as  being 


«(x)  -  (t)  =  21  (T)  +  v(T) 


(4.  17) 

(4.  18) 
(4.  19) 
(4.  20) 

(4.  21) 

(4.  22) 

(4.  23) 
(4.  24) 

(4.  25) 

(4.  26) 
of 

(4.  27) 
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The  above  linearization  procedure  yields  the  following  results.  The 
linearization  of  the  nonlinear  system  with  x  as  the  state  vector  leads  to  a 
linear  and  time -varying  system,  L^,  with  state  vector  £  and  output 
vector  r\  related  by 


L 


1 


F*(t)£ 

H*(t)  5 


(4.  28) 
(4.  29) 


On  the  other  hand,  the  linearization  of  the  nonlinear  system  with  z  as  the 
state  vector  leads  to  a  linear  and  time -varying  system,  L^,  with  state 
vector  J  and  output  vector  p  related  by 


1  =  F2(t)j 

=  H*(t)J 


(4.  30) 
(4.  31) 


We  shall  now  proceed  to  derive  the  relations  which  relate  the  matrices 
F I  (t) ,  H*(t),  G*(t),  F*(t),  and  H^(t)  in  Eqs.  (4.  28)  through  (4.  3  1 ).  To  do 
this,  we  start  from  Eq.  (4.4),  namely, 


z'  (t)  =  j;(x'‘'(t)  ) 


(4.  32) 


Therefore, 


z"(t)+  J(t)  =  (t)  +  J(t) 


=  _g(x '  (t)  )  + 


Ox 


x‘  (t) 


?(t)  +  o2(§) 


(4. 33) 
(4.  34) 


Let  (see  Eq.  (3.  7)  ) 


G  (t)  2  G  ( x  ( t)  )  g 


X  (t) 


(4.35) 


Then,  from  Eqs  (4.35),  (4.34),  and  (4.32)  we  deduce  the  relation 


J(t)  =  GV(t)I(t)  +  o2(f) 


(4.36) 
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Neglecting  the  higher  order  terms  c>2(jj  we  ^ave  f°r  sma-H  j|  _|(t) 


J(t)  =  G  (t)_£  (t) 

Differentiating  with  respect  to  time  we  have 
J(t)  =  G*(t)?(t)  +  G*(t)j(t) 

Substituting  Eqs.  (4.  3  0)  and  (4.  28)  into  Eq.  (4.  3  8)  we  obtain 


F2(t)J(t)  = 


G  (t)  +  G  (t)  F  j  (t) 


§(t) 


(4. 37) 


(4.38) 


(4.  39) 


Substituting  Eq.  (4.  37)  into  Eq.  (4.  39)  and  noting  that  the  ensuing  relation 
must  hold  for  all  §  (t)  we  deduce  that 


F.2  (t)  G  (t)  =  G  (t)  +  G  (t)Fj  (t) 


or 


,,-1 


*-l 


F2(t)  =  G  (t)  G  (t)  +  G  (t)  Fj  (t)  G  (t) 


(4.  40) 


(4.41) 


Clearly,  Eq.  (4.41)  links  the  matrices  F2  (t)  and  F  (t)  appearing  in  Eqs. 
(4.  28)  and  (4.  30). 

5|C  5|< 

Next,  we  proceed  to  relate  the  matrices  Hj  (t)  and  (t).  From 
Eqs.  (4.29),  (4.31),  and  (4.  3  7)  we  conclude  that 


Let 


H*(t)  =  H^tJG^'V) 


❖  A  *-l 

r  (t)  =  g  (t) 


Then  Eqs.  (4.  41)  and  (4.  42)  reduce  to 

F*(t)  =  G*(t)£*(t)  +  G*(t)F*(t)r*(t) 

H2(t)  =  H*(t)r*(t) 


(4.42) 

(4.43) 

(4.  44) 
(4.  45) 


We  remark  that  the  above  equations  hold  if  and  only  if  the  various  higher 
order  terms  are  neglected. 


5.  THE  KALMAN  FILTERS  AND  THE  COVARIANCE  EQUATIONS 


Let  us  consider  the  linear  system  Lj  defined  by  Eqs.  (4.  28)  and  (4.  29). 
The  observed  signal  is  (see  Eq.  (4.27)  ) 


w(t)  =  ri(t)+v(t)  (5.1) 

We  can  construct  a  Kalman-Bucy  filter  (see  Ref.  )  which  will  yield  an 
estimate  §  (t)  of  £(t).  The  filter  is  specified  by  the  equation 

F*(t)  -  R”  1(t)  £(t) 

+  ;  £(to)  =  0  (5.2) 

A 

where  £_(t)  is  the  covariance  matrix  of  the  error  vector  £ (t)  -  £(t)  and 
it  satisfies  the  Riccati  matrix  differential  equation 


d_ 

dt 


I'. 


Fj  (t)£ 


-§  +-§-l 


F,  (t) 


£1]  (*)£§ 


(5.3) 


Similarly  consider  the  linear  system  L^  defined  by  Eqs.  (4.  3  0)  and 
(4.31)  with  the  observed  signal  given  by  Eq.  (5.  1).  Again  we  can  construct 
a  Kalman-Bucy  filter  which  generates  an  estimate  J(t)  of  J(t).  The  filter 
is  specified  by 


=  ^F*(t)  -  £j(t)H*V)R'^t)H*(t) 


J(t) 


+  £J(t)H*'(t)R'1(t)£(t)  ;  J(tQ)  =  0 


(5.4) 


where  j^j(t)  is  the  covariance  matrix  of  the  error  vector  J^(t)  -  J(t)  and 
it  satisfies  the  Riccati  equation 


£  F 
-J  -2 


(t) 


-  £  XX 

— J  —2 


H*'(t)R-1(t)  H 


-2 


(t)£j 


(5.  5) 


We  shall  now  show  that  the  covariance  matrices 
related  by  the  relation 


£^(t)  and  £  (t)  are 
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(5.6) 


E,(t)  =  r:,;(t)^j(t)r;'c  (t) 

To  do  this  wc  proceed  as  follows.  Let 


xu)  § 


Then, 


~ ;  X(t)  =  r(t)^(t)T'(t) 

+  r(t)^j(t)_r'  (t) 

Substituting  Eq.  (5.  5)  into  Eq.  (5.  8)  we  obtain 

x(t)  =  + 

I  (t)  — J(t/  — Z 


(5.  7) 


(5.  8) 


-  H2* '(t)  R'  *(t)  H* (t)  Ij(t) 


(5.9) 


Next  we  substitute  Eqs.  (4.44)  and  (4.45)  into  Eq.  (5.  9);  we  use  the  relation 

vl/  J/ 

that  (since  from  Eq.  (4.  43),  £  (t)  G  (t)  =  I_  ) 

£*(t)  =  -£*(*)  G*(t)£*(t)  (5.10) 

to  obtain 

x(t)  =  -i*(t)6*(t)  L*(t)fj(t)  r'(t) 

-  r;(t)£J(t)r,(t)Gc'(t)r,(t) 

+  r*(t)G*(t)  r*(t)Ej(t) 

+  I*(t)G*(t)F*(t)I*(t)£J(t)r*,(t) 
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Substituting  Eq.  (5.  7)  into  Eq.  (5.  11)  we  obtain 


X(t)  =  +  X(t)  F*\t)G*'(t)r*'(t) 

-  X(t)H*’(t)R'1(t)H*(t)X(t)  (5.12) 

Recalling  that  I  (t)  G  (t)  =  I_  we  see  that  Eq.  (5.  12)  reduces  further  to 

X(t)  =  F j  (t)  X(t)  +  X(t)  F  j  ' (t)  -  X(t)H*V)R_1(t)H*(t)X(t)  (5.  13) 

Let  us  now  compare  Eqs.  (5.  13)  and  (5.  3).  We  conclude  that  X(t)  and 
£-(t)  satisfy  the  same  matrix  differential  equations.  Furthermore,  from 
Eqs  (5.7)  and  (4.37),  X(to)  =  Z*(t0)^j(t0)£*'(t0) 

X(tD)  =  I*(t0)£J(t0)l’>,,‘0) 

=  iX> E  {kvj'1 (‘o>}  r’(to» 

=  f  (‘0)G*(‘o>  E  S*'(t0)r*'(to) 

=  E(i(t0)i'(t0)J  4  %t0) 


(5.  14) 


so  that  X  and  have  the  same  initial  conditions.  Therefore,  by  the 

uniqueness  theorem  of  solutions  to  differential  equations  we  conclude  that 


X(t)  =  £q*(t)  for  all  t 


(5.  15) 


and,  hence,  from  Eq.  (5.  7),  that 


Up  to  now  wo  have  shown  that  the  two  error  covariances  matrices 
^r(t)  and  a  re  related  by  the  linear  equation  (5.  16);  this  can  also 

be  written  as 

ilj(t)  =  G  (t)  £-(t)  G  '  '(t)  (5.]  6) 

A  r 

The  next  step  is  to  compare  the  estimates  j_(t)  and  ^J(t)  of  the  two 
Kalman  filters  (5.  2)  and  (5.  4).  We  shall  show  that 

J(t)  =  G*(t)|(t)  (5.17) 

To  do  this,  define  the  vector 

S(t)  =  r(t)J(t)  (5.18) 


Note  that 


q(t  )  =  1  (t  )  J  (t  )  =  0  =  ?  (t  ) 
_  o  —  o  —  o  —  —  o 


(5.  19) 


We  shall  show  that  q(t)  and  §  (t)  satisfy  the  same  differential  equation; 
also,  by  (5.  19)  they  have  identical  initial  conditions;  this  would  imply  that 

A 

q(t)  =  _£(t)  and,  therefore,  Eq.  (5.  17)  would  follow. 

From  Eq.  (5.  18)  we  have 


d 

dt  2J 


q(t)  =  f(t)l(t)  +  r*(t)  ~J(t) 


=  -T(t)G,!(t)_rv(t)  J(t)  +  rv(t)  ^J(t) 

=  _  r*(t)G*(t)  q(t) 

+  [ £*(t)  F  * (t)  -  £*(t)  L  j(t)  H* ’(t)  R_I(t)  H*(t)]  J(t) 
+  £>:<(t)  2  j(t)H* ’(t)  R_1(t)£(t) 


(5.  20) 

(5.  21) 
(5.  22) 


(5.  23) 


Substituting  Eqs.  (4.44),  (4.45),  and  (5.  16)  into  Eq.  (5.23)  we  obtain 
after  some  algebraic  manipulations 
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(5.  24) 


£  q(t)  =  [F*(t)  -^t)H*'(t)R“1(t)H*(t)]  q(t) 

A 

Comparing  Eqs.  (5.24)  and  (5.  2)  we  deduce  that  J^(t)  and  q  ( t )  satisfy 
the  same  differential  equation.  Thus,  on  the  basis  of  our  previous  argument 
we  conclude  that 

J(t)  =  G*(t)|_(t)  (5.25) 

The  implications  of  this  result  will  be  discussed  in  the  following 
section. 
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6.  Conclusions 


Let  us  recapitulate  the  development  and  the  results  up  to  now. 

We  started  with  two  different  state  variable  representations  of  the  same 
system  (see  Eqs.  (2.1),  (2.2)  and  (2.5),  (2.6)).  We  linearized  the  system 
about  the  same  nominal  trajectory.  In  this  manner  we  derived  two  linear 
systems,  describing  the  equations  of  the  perturbation  vectors,  L^  and  L^  (see 
Eqs.  (4.28)  through  (4.31)).  These  systems  were  described  by 

i(t)  =  F*(t)  £  (t) 

J(t)  -  A*(t)  J(t) 

We  showed  that  (see  Eq.  (4.37)) 

J(t)  =  G*(t)  §  (t)  (6.2) 

to  within  first  order  terms  where  G  (t)  is  a  computable  matrix  from  the 
problem  data  (see  Eq.  (4.35)). 

Next  we  showed  how  to  generate  the  estimate  ^(t)  of  (t)  and  the  estimate 

<v 

J(t)  of  J(t)  using  the  Kakman-Bucy  Filter.  After  many  manipulations  we 

/v  A 

showed  that  the  estimates  (t)  and  J^t)  were  related  by  (see  Eq.  (5.  25)) 

,  /V 

J(t)  =  C  (t)  §  (t)  (6.3) 

In  other  words,  the  estimates  are  related  by  the  same  linear  transformation, 
G  (t),  as  the  perturbation  vectors.  Since  G  (t)  is  a  deterministic  matrix  we 
conclude  that  we  can  construct  a  Kalman  Filter  based  upon  the  linearized 
state  equations  in  any  one  coordinate  system  to  obtain  the  estimate,  _x  (t),  of 
the  state  x(t)  in  that  coordinate  system  by 

/\ 

x(t)  =  x"(t)  +  S  (t)  (6.4) 


(6.1) 
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where  x  (t)  is  the  deterministic  nominal  trajectory.  If  we  chose  a  differ¬ 
ent  coordinate  system  in  which  the  state  is  z(t)  and  given  that  _z  and  x 
are  related  by 


z  =  g  (x)  or  x  =  Y  (z) 


(6.  5) 


Then  if  we  let 

*  A  ! 

5  (t)  =  l"3x"/x*(t) 

then  we  can  construct  the  estimate  _z(t)  by 
z(t)  =  g(x'"(t)  )  +  GV(t)i(t) 


(6.  6) 


(6.7) 


using  the  output  of  the  previous  Kalman  filter.  If  higher  order  terms  in 
the  linearization  scheme  are  negligible,  the  estimate  ^(t)  is  optimal. 

The  above  discussion  does  not  mean  that  there  need  not  exist  a  co¬ 
ordinate  system  in  which  the  linearization  is  more  accurate  than  in  another 
coordinate  system.  However,  it  is  not  known  how  to  find  such  a  coordinate 
system. 
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